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ABSTRACT 

Sensitive and high angular resolution (~ 0.4") S02[222,2o 22i 21] and SiO[5— >4] 
line and 1.3 and 7 mm continuum observations made with the Submillimeter Array 
(SMA) and the Very Large Array (VLA) towards the young massive cluster W5 1 IRS2 
are presented. We report the presence of a large (of about 3000 AU) and massive (40 
Mq) dusty circumstellar disk and a hot gas molecular ring around a high-mass proto- 
star or a compact small stellar system associated with W51 North. The simultaneous 
observations of the silicon monoxide molecule, an outflow gas tracer, further revealed 
a massive (200 M0) and coUimated (~ 14°) outflow nearly perpendicular to the dusty 
and molecular structures suggesting thus the presence of a single very massive proto- 
star with a bolometric luminosity of more than 10^ L©. A molecular hybrid LTE model 
of a Keplerian and infalling ring with an inner cavity and a central stellar mass of more 
than 60 M0 agrees well with the S02[222,2o — > 22i 21] line observations. Finally, these 
results suggest that mechanisms, such as mergers of low- and intermediate- mass stars, 
might be not necessary for forming very massive stars. 

Subject headings: stars: pre-main sequence - ISM: individual: (W51 IRS2) - ISM: 
Molecules - radio continuum 
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1. Introduction 



In recent years a sm all group of candidate accreting disks in high-mass protostars has been 



reported in the literature (|Cesaroni et al 



120061 and references therein), but with luminosities 
typical of B-type main-sequence stars, that is, stars with masses < 25 Mq. Around more luminous 
and massive objects (presumed 0-type stars) there have been no clear evidence for accretion 
disks, only gravitational unstable and large rotating molecular structures have been found. These 



molecular structures i^'toroids") are infalling and thus accreting fresh gas rnateria 



cluster of young massi ve protostars ( SoUins & Ho 



single 40 M© protostar (|Sandell et al. 



2003 



2005 



Beltran et al. 



Beuther & Walsl 



3 20! 



2005L 



to a central 



20061) or maybe a 



The sizes and masses of the 



rotating toroids are about 2-3 x 10"^ AU and 100-400 Mq, respectively. 



2002 



ocated at 5-8 kpc a way in the Sagittarius spiral arm (iGenzel et al.lll98ll : llmai et al. 



Schneps et al. 



19811) and with a total bolometric luminosity of about 3 x 10^ Lq, the 



W51-IRS2 region is one of the most luminous massive star forming regions in our Galaxy 



(lErickson & Tokunaga 



1980). 



Recently however it has been es timated a distance of 2.3 + 0.3 kpc to W51-IRS2 using 



spectroscopic parallaxes of OB stars (|Figueredo et al 



20081 ). But, in other hand, a lower limit of 



5 kpc by triangulation using Very Long Bas eline Array observations of methanol masers located 



in W51 IRS2 was found bv ,Xu et al 



statistical H2O maser parallax (IGenzel et al. 



(I2OO8I). a ve r y similar value to those found earlier by the 



1981 



Schneps et al. 



1981 



Imai et al 



20021) . In this 



study, we adopt a distance of 6 kpc for the cluster. 

This region might contain about 30 0-type zero-age main-sequence star s with excess 



Barbosa et al 



emissi on at the infrared and (sub)millimeter wavelengths, see for an example 
(l2008b . In the W51 IRS2 cluster there are two highly obscured massive sources called W51 
North and W51 d2 that seem to be the yo ungest objects i n the cluster and tha t exhibit strong 



emission from many complex molecules (IHo et al 



1983 



Rudolph et al 



19901 : 



Zhang & Ho 
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19971 : IZhang et alJI 19981 : ISollins et alJl2004l : IZapata et al.ll2008bh . In particular, toward the W5 1 
North object very strong maser emission from vario us species e.^. , hydroxy! (OH), silico n 



monoxide (SiO), a nd water (H2O) has been detected dSchneps et al. 



1981 



Gaume & Mutel 



Morita et al. 



19871 : 



1992|) . This type of emission is associated with the formation of the high-mass 



stars. Observations of the proper motions of the H2O masers, in addition, revealed the presence 



of a remarkab 



dSchneps et al 



e cornpact (~ 1-2") southeast- northwest ou tflow emanating from this protostar 



1981 



Eisner et al. 



2OO2I : 



Imai et al 



20021) • Furthermore, very high spatial 



resolution VLA and VLBA radio observations revealed that the SiO maser emission traces the 



innermost parts of the outflow ejected from W51 North (lEisner et al. 



20021) • Finally, observations 



of the cyanogen (CN) molecule showed th at the molecu 



with a mass accretion rate ~ 10 ^ Mq yr ^ (IZapata et al. 



ar gas is falling into this central protostar 



2008b ). 



2. Observations 



In order to study the very young massive object W51 North, we carried out observations of 
the SO2 [222,20 ~^ 22i 21] (sulfur dioxide) and SiO[5^4] (silicon monoxide) molecular lines and 
1 .3 mm continuum with the Submillimeter Array. We also obtained 7 mm continuum data with 
the Very Large Array. Both observations had high angular resolution (~ 0.4"). 



The SO2 molecule i s a good high-density gas tracer (see 



Leurini et al. 



Jimenez-Serra et al. 



20071 : 



20071), while the SiO, is an excellent tracer for molecular outflowsl 
The SMA SiO[5^ 4], SO2[222,20 22i 21] and 1.3 mm observations were made simultaneously 
on 2008 January 17, and the VLA 7 mm observations were made on 2008 April 22. At that time. 



the SMA was in its very extended configuration, while the VLA was in its C configuration. The 



' The production of the silicon mo noxide (SiO) molecule is mainly attributed to the destruction 
of the dust particles in strong shocks (jSchilke et al.ll 19971 : Ivan Dishoeck & Blakelll998|) . 
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phase reference center of both VLA and SMA observations was R.A.=19'' 23™ 40'.'055, decl.= 14° 
31'5'/59(J2000.0). 

For the 7 mm observations, we integrated on-source for a total of approximately 5 hr using 
the fast- switching mode with a cycle of 120 s. The observations were made using the continuum 
mode, with a total bandwidth of 100 MHz. The central frequency observed was 43.34 GHz. The 
absolute amplitude calibrator was 1331+305 (with an adopted flux density of 1.45 Jy) and the 
phase calibrator was 1924+156 (with a bootstrapped flux density of 0.65 + 0.01 Jy). 

For the 1 mm observations, the zenith opacity (t23oghz), measured with the NRAO tipping 
radiometer located at the Caltech Submillimeter Observatory, was ~ 0.15, indicating excellent 
weather conditions (for this frequency) during the observations. Observations of Uranus provided 
the absolute scale for the flux density calibration. Phase and amplitude calibrators were the 
quasars 1925+211 and 2035+109, with measured flux densities of 0.7 + 0.1 and 0.5 + 0.1 Jy, 
respectively. 

The SO2[222,20 ^ 22i 21] and SiO[5^4] transitions were detected in the lower side band 
(LSB) of the SMA at a frequency of 216.643 and 217.104 GHz, respectively. The fuU bandwidth 
of the SMA correlator is 4 GHz (2 GHz in each band). The SMA digital correlator was configured 
in 24 spectral windows ("chunks") of 104 MHz each, with 256 channels distributed over each 
spectral window, providing a resolution of 0.40 MHz (0.58 km s"^) per channel. Fur ther technical 



descriptions of the SMA and its calibration schemes can be found in , Ho et al. 



(120041) 



The data were edited, calibrated and imaged using the programs MIR, AIPS, MIRIAD 
and KARMA. In both bands, we used the ROBUST parameter set to to obtain an optimal 
compromise between sensitivity and angular resolution. The 7mm data were self-calibrated in 
phase. For the 7 mm observations, the continuum image rms noise was 0.9 mJy beam^ at an 
angular resolution of 0!'46 x 0"43 with a P.A. = -53.6°. For the 1 mm observations, the continuum 
image rms noise was 6 mJy beam"^ at an angular resolution of 0"58 x 0"43 with a P.A. = 54.8°. 
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The line image rms noise was about 80 mJy beam 



3. Results and Discussion 



3.1. A massive and large dusty disk 

At a wavelength of 7 mm (see Figure 1), we detected strong free-free continuum emission 
arising from the HII regions associated with the neighboring young massive ZAMS stars located 



in the W5 1 IR S2 cluster. This emission has already been mapped at centimeter wave 



many authors (iGaume & Mutel 



1987; 



Eisner et al. 



radio images presented by ,Lacy et al. 



20021 : 



Mehringei 



1994; 



Lacv et al. 



2007|) were obtained from , Mehringei] (|l994|) . However, 



engt 



IS by 



2007). The 



in our map, we detected for the first time a faint and compact millimeter source associated with 
the obscured high-mass protostar located in W51 North. This 7 mm continuum source is the 
counter part of the d usty sources associate d with large envelopes, and reported at 2 and 1.3 mm by 



Zhang et al. 



(1998); 



Zapata et al. 



(l2008bh . respectively. 



In Figure 2, we have constructed the spectral energy distribution (SED) for this millimeter 
source from the centimeter to millimeter wavelengths, combining only data presented here. The 
angular resolution of two different observations are quite similar, of about 0.4". From these data, 
we estimated a spectral index of or = 2.8. This suggest that the emission at these wavelengths 
likely arises from a dusty disk. A larger envelope has been already mapped at scales lO'^ AU 
as mentioned above. The hypothesis of a flattened circumstellar disk is supported by the fact 
that the mm source is resolved at wavelength 7 mm (see Table 1) and show a modest east- west 
elongation, nearly perpendicular to the orientation of the molecular outflow that will discussed on 
the following sections. At 1.3 mm the orientation is not so well determined. 



We noted from Figure 16 of 



Mehringeii (119941) there is no strong emission from the 



recombination lines H92Qr and He92Q' toward the mm source, both lines are associated with the 
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compact HII regions W51d and/or dl. The non detection of recombination lines emission from 
the mm source is consistent with our interpretation. Taking the values of the deconvolved sizes 
and the flux densities for the continuum source at 7 mm and 1.3 mm, we estimated a brightness 
temperature 80 + 10 K., a very low temperature to be associated with a H II region. This 
temperature is more likely associated with thermal dust emission. 

Possibly, the high mass accretion rates of the infalling material associated with W51 North 
(~ 10"^ Mq yr"^) could be likely quenching or trapping the development of an ultra-compact 



HII region, so that the free-free emission from the io nized materia^ 



wavelengths, see for a reference of these phenomena: 



Osorio et al. 



is undetectab 



mm; 



e at centimeter 



Ketol (120031) . Another 



possibly is we have a quite young massive protostar that has not developed an HII region yet. 

In Figure 1, we also show the morphology of the continuum source at 1.3 mm. At both 
wavelengths, the deconvolved sizes are similar, with a size on the order 3500 x 1500 AU, see 
Table 1. Certainly, a very large disk compared with those observed in low-mass stars (~ 100 



AU), but with similar sizes to those observed in early B 



- ^— /? ■■ 

Cesaroni et al. 


1999; 


2007 


). 



Shepherd et al. 



2001 



type protostars ( 



Schrever et al. 



20061 : 



Patel et al. 



000-2000 AU, e.g. 



2005 



Rodriguez et al. 



Finally, 



dZhang et al. 



assum ing that the emission is optically thin, a dust temperature value of 200 K 



19981) . a grain emissivity spectral index j3=\ (see Figure 2), a value of the dust 
absorption coefficient that goes as [ ^.Jl^i ] = ^-^I- ioooghz ^^ ^^'^ ^ distance to W51 North of 6 kpc, 
we estimated a gas mass of the dusty disk of 40 Mq. This value is cons istent with the mass 



(100 Mq) obtained f or the large dusty structures found in W5 1 North by 



Zapata et al. 



(2008b) 



Zhang etal.l(ll998h : 
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3.2. A coUimated molecular outflow 



The SiO[5— >4] zero moment map shows the presence of a compact and collimated (~ 14°) 
north- south bipolar molecular outflow. The redshifted radial velocities are from +60 to +95 km 
s'^ and the blueshifted ones are from +20 to +58 km s"^; the systemic LSR radial velocity of 
the ambient molecular cloud is about at 59 km s^^ The receding radial velocities (redshifted) 
of the outflow are located toward the north, while the approaching radial velocities (blueshifted) 
are toward the southwest. The outflow emanates from the mm source imaged at 1.3 and 7 mm 
wavelengths (see Figure 1). 

In Figure 1, we have also overlaid the SiO[5^4] molecular emission fro m the bipolar outflow 



with the positions of the blue- and red-shifted water maser spots reported by 



Eisner et al. 



(2002) 



It is evident how the water maser spots trace the innermost parts of the SiO[5^4] molecular 



bipolar outf 



Zapata et al. 



ow as o bserved on many other molecular outflows and proposed for this case by 



(2008b). 



In the other hand, the SiO masers reported also in 



very luminous high-mass star forming regions (|Zapata et al 



Eisner et al. 



( 20021) and associated with 



2008ah . might be tracing even 



more dense parts of the outflow at very small scales, as has been observe d for the massiv e 



protostar Source I located in the closest massive star forming region, Orion (|Menten & Reid 



Reid et al 



1995 



20071). We have marked the position of the SiO maser source in our Figure 1. The RA. 



of the SiO[5^4] molecular outflow is 150° + 20°, almost perpendicular to the orientation of the 
dusty source, as mentioned before. The large radial velocities of the outflow and the morphology 
of the SiO[5^4] molecular emission suggest that it is nearly perpendicular to the plane of the sky. 
Thi s is also suggeste d by the location of water maser radial redshifted and blueshifted velocities. 



see 



Imai et al. 



(l2002h : 



Eisner et al. 



(2002). 



Assuming that we are in LTE and the molecular emission is optically thin, we can estimate 
the mass of the outflow for the SiO molecule in the transition AJ= 5^4 from the equation: 
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Mh^ = 1.75 X 10 



Xs^ov^ (l-exp[^]) 



(1) 



where j^, is the fractional abundance between the silicon monoxide and the molecular hydrogen, 
that for this case we took a value of 10^, D is distance to W51 North, AQ. is solid angle, v the 
rest frequency of the line, Sy the flux density and the excitation temperature equal 50 K. 
We estimated a mass for the outflow of 200 Mq. It is important to mention that the fractional 
abundance between the silicon monoxide and the molecular hydrogen varies in every star forming 
regions (SFRs) and thus the mass of the outflow can be over/underestimated. The value for the 



fractio n al abundance assumed here seems to be c onsistent for a few SFRs, see .Ziurys & Friberg 



ll987t); 



Mikami et al. 



(1992); 



Zhang et al. 



( 1995 ). 



For a mechanical force of F m=200 Mq 30 k m s / 1800 yr =3.3 Mq km s yr and 



from the correlation presented in 



Wu et al. 



( 20041) for the outflow mechanical force versus the 



bolometric luminosity of the exciting source, we very roughly estimated a luminosity for the 
central powering source of more 10^ Lq which corresponds to a massive early 0-type protostar. 
This spectral type for the central star is in good agreement with that obtained from the dynamical 
considerations. 



3.3. A hot molecular gas ring 

A large and flattened molecular ring (with an inner cavity of about 3000 AU) surrounding the 
dusty disk was traced by the emission of the molecule S02[222,2o 22i 21] (see Figure 1). 

The ring has a deconvolved radius of order 6000 AU with an orientation (RA.= 22°) almost 
perpendicular to that of the molecular outflow and similar to the presumed dusty disk. In Figure 
3, we show the position-velocity diagram of the molecular S02[222,2o — > 22i 21] emission from the 
ring structure. From this figure, one can see a ring structure with a modest velocity gradient across 
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the ring of ~ 5 km s"^ over 2.0", produced by rotation and infall. The redshifted radial velocities 
located to the east, while the blueshifted ones to the west. The small velocity gradient is likely due 
to the small inclination angle of ring with respect to plane of sky. Similar rings velocity structures 
were also found traced by the emission of other molecules species, e.g HC3N, H2CO, NH2CHO 
and CH3OH. 



Very recently, however, IZapata et al.l (l2008b|) reported a large structure or envelope with 
a size of 4 x 10"* AU and with a P. A. = 90° traced by the C2H3CN molecule and centered in 
W51 North. This molecular structure seems rotate in an opposite direction (with the blueshifted 
velocities to the east while the redshifted ones to the east) to the mo lecular compact ring traced 



by the SO2[222,20 — > 22i_2i] . As discussed in 



Zapata et al. 



(|2008bh . other molecules such as 



HCOOCH3 and CH3OH at same scales associated with W51 North were very much contaminated 
by the emission from the hot molecular core associated with W51d2 do not allowing to confirm 
the velocity shift. We thus think that the C2H3CN molecular core in W51 North was also 
contaminated with emission coming from W51d2 given therefore the impression of rotation. 

We are more confident of the velocity gradient across the molecular ring are correct because 
we have resolved the molecular source and other molecules show similar velocity structures. 

In Figure 4 and 5, we additionally present the velocity channel maps and the position- velocity 
diagram of the S02[222,2o — ^ 22i 21] emission from the ring, overlaid with a simple LTE model of 
a large infalling Keplerian ring with an inner cavity with sizes (of the ring and the inner cavity) 
similar to those observed. 

We assumed in this model that the thermal molecular line emission of the ring is in local 



thermodynamic equilibrium, and took from the lite rature some 



mow values as the excitation 



Zhang et al. 



1998 



Zapata et al. 



2008b). The 



temperature, the density and the infall velocity (see 
best fit was found recurrently until we obtained similar structures in our model to those imaged 
(see. Table 1 and Figures 3 and 4). The model fits the observations well. 
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This model, in addition, restricted the position of the ring on the sky. It is almost in face-on, 
with an inclination angle with respect to plane of the sky less than 30° and a P. A. equal to 22°. It 
is important to mention that if the inclination angle is smaller, the dynamical mass will increase 
by much more, with a reason inversely proportional to the sin^(/), where i is the inclination angle. 
Thus, the central object could have a stellar mass of more than 60 M©. 

It is interesting to note that the mass of the disk (40 Mq) might be on the order of the stellar 
mass (> 60 Mq). This has seen before for the massive young disks associated with the high-mass 



protostars located in the star forming region: NGC6334N(I) (Rodriguez et al. 



2007h 



The nature of large inner cavity in the molecular ring may be due to opacity effects, 
photodissociation of the molecules, or clearing of the inner disk either from a multiple system of 
compact circumstellar disks already formed in the center of the ring or by the outflow itself. If the 
nature of the ring is because of opacity eff'ects or photodissociation of the molecules very close to 
the central object suggesting that both the molecular ring and circumstellar disk could be part of 
a same extended and flattened structure. Observations with a much better angular resolution of 
high-density tracers and the continuum at millimeter wavelengths may determinate the nature of 
the inner cavity in the ring. 

The strong feature observed toward the east of the ring might be originated from the 
interaction of the outflow with the ring itself. However, more observations at high angular 
resolution with different molecular probes sensitive to different cloud properties will be quite 
important for studying the higher excitation material within the ring. 

In Figure 5 we show artist's conception of the molecular ring, the bipolar outflow and the 
central massive disk found in the high-mass pro to star W51 North. 
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4. Summary 

W51 North presents a promising laboratory for future studies on formation of very high-mass 
stars. We have observed this young massive protostar at spatial scales of 0.4" using the SMA and 
VLA and found the following: 

• A massive (40 Mq) and large (~ 3000 AU) dusty circumstellar disk around a single 
high-mass protostar or a compact small stellar system. 

• A hot and rotating molecular ring with a inner cavity of about 3000 AU and a size of about 
9,000 AU. 

• A highly coUimated and massive (200 Mq) bipolar outflow that emanates with an orientation 
almost perpendicular to the circumstellar disk and the ring. This coUimated outflow suggests 
the presence of a single massive protostar in the center of the circumstellar disk. 

• Our data reveled that H2O and SiO masers associated with this highly embedded protostar 
are tracing the innermost parts of the SiO [5 ^4] thermal bipolar outflow. 

• A molecular hybrid UTE model of a Keplerian infalling ring with an inner cavity (~ 3000 
AU) and a stellar mass of more than 60 M© agree well with the SO2[222,20 22i,2i] line 
observations. 

These results suggest that mechanisms, such as mergers of low- and intermediate- mass stars, 
might be not necessary for forming very massive stars. 

We thank all the SMA and NRAO staff" members for making these observations possible. 
The Submillimeter Array (SMA) is a joint project between the Smithsonian Astrophysical 
Observatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded 
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Table 1. Physical Parameters of the Molecular Ring and Circumstellar Disk 



Position 



Wave./Mol. ff2000 ^2000 Flux Density Deconvolved Size P.A. Gas Mass Dyn. Mass 

[/!)»,v| |o///] |mJy] larcscc] IDcgrccl |\1 ] |Mv| 

7 mm 19 23 40.057 14 31 05.67 17 ± 3 0.58±0.02 x 0.27±0.02 70±6 

1.3 mm 19 23 40.045 14 3105.48 1620 ± 160 0.72±0.03 x 0.57±0.03 90±40 40 

SO2 19 23 40.074 14 31 05.44 7x10* ± 500^ 1.8±0.1 x 1.2±0.1 72±2 - > 100^ 

Note. — 1: The units are mjy Beam"^ km s'^. 2: The dynamical mass was obtained from our model. 

Table 2. Parameters for the molecular gas ring LTE model 



Name Parameter Value 



Systematic Velocity 
Orientation 
Inclination 

Density 
Iimer Radius 
Outer Radius 
Temperature of the central object 
Dust Exponent 

Distance 
Reference radius 
Power law index density 
Kinetic temperature 
at referenced radius 
Excitation temperature 
at referenced radius 
Power law index Tfo„ 
Power law index 
Scale Height of disk 
Dynamical mass 
infall velocity 



Ylsr 60.0 km S-' 

PA 22° ± 10° 

/ < 30° 

p 10* cm-3 

R; 3000 AU ± 200 AU 

Ro 6000 AU ± 200 AU 

T 500 K 

/? 1.0 

D 6000 pc 

r 6000 AU 

a 2.75 

Tfa„ 200 K 

Texc 200 K 

y 0.6 

6 0.6 

H 60 AU ± 20 AU 

M > 100 Mq 

V,„/ 4kms"' 
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by the Smithsonian Institution and the Academia Sinica. NRAO is operated by Associated 
University, Inc., under contract with the National Science Foundation. 

Facilities: Submillimeter Array (SMA), Very Large Array (VLA). 
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7000 AU. ' v'- W51d2^'^^^ 

. - * W51 North 



40,00^ 

Right Ascension (J2000) 

14°31'06.5' 





19''a3™40.20" 



40.10" 40.00" 
Right Ascension (J2000) 



9^a3'°40. 12040, 100^ 40.080^ 40.060" 40.040" 40.020" 
Right Ascension (J3000) 



Fig. 1 . Upper Panel: Overlay of an archival VLT H-band image (grey scale), a VLA 7 mm continuum image (grey contours) and a SMA 

SiO[5— >4] moment zero line image (blue and red contours) of the young massive cluster W51 IRS2. Blue and red contours correspond respectively 
to the blue- and red-shifted gas that emanates from the obscured protostar W51 North. The 7 mm contours are from 10% to 90% with steps of 3% 
of the peak of the continuum emission, the peak is 100 mjy beam"'. The synthesized beam of the 7 mm contour image is 0.46" X 0.43" with a P. A. 
of -53.6° and it is shown in the upper right comer of the box. The SiO[5— »4] line contours are from 20% to 90% with steps of 10% of the peak of 
the line emission, the peak is 5 Jy beam"' km s"' for the red-shifted emission and 8 Jy beam"' km s"' for the blue-shifted. The integrated velocity 
range for the blueshifted gas is from +20 to +58 km s"' and that for the redshifted gas is from +60 to +95 km s"' . Left Lower Panel: A zoom into 
the W51 North zone. Overlay of the SO2[222,20 ~^ 22i 21] moment zero line emission (color image), the 7 mm continuum emission from the disk 
(white contours) and the SiO[5— >4] moment zero line emission. The integrated velocity range for the SO2[222,20 ~> 22i 21] is from +50 to +70 km 
s"' . The synthesized beam of the SiO[5— >4] contour image is 0.58" X 0.43" with a P.A. of 54.8° and it is shown in the bottom right corner of the 
box. Right Lower Panel: A zoom into the dusty disk and the outflow. The grey scale image shows the 1.3 mm line-free continuum emission overlaid 
with the 7 mm continuum emission (white contours) and SiO[5— >4] line emission (blue /red contours). The blue and red dots mark the position of 
the blue- and red-shifted water maser spots, respectively, reported by|j 



Eisner et al. 



20021) . The synthesized beam of the 1.3 mm continuum image is 



0.50" X 0.35" with a P.A. of 54.8° and it is shown in the bottom right comer of the box. For the 1.3 mm image, we used the ROBUST parameter of 
the INVERT task of MIRIAD set to -2, w hich corresponds to uniform weighting to achieve the maximum angular resolution. The black dot marks 
the position of SiO maser source given bv Eisner et al 



2OO2I) , (a, iJ)j2000= (19* 23'" 40'.'055. 14° 31' 5"59). 
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Fig. 2. SED for the source W51 North combining 7 and 1.3 mm continuum data. The respective error bars were smaller than the squares 

and are not presented. The line is a least-squares power-law fit (of the form Sy <x v") to the spectrum. 
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Fig. 3. Velocity-channel image of the SO2[222,20 —> 22i_2i] thermal emission from the ring (contours) overlaid with our synthetic LTE 

model (greysacle). The spectral resolution was smoothed to velocity bins of 1 km s"^ The central velocity is indicated (in km s"') in the top 
rigth-hand comer of each panel. The systemic velocity of the ambient molecular cloud is about 60.0 km s"^ The contours are 10% to 90% with a 
step of 10% of the peak flux. The scale bar indicates the peak flux of the molecular emission in T;,. 
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Angular Offset 

Fig. 4. Position-velocity diagram of the molecular ring computed at P.A. = 22° (colors) overlaid with the position-velocity diagram from 

our model (contours). The contours are from 10% to 90% with steps of 10% of the peak of the line emission of our LTE model. The units of the 
vertical axis are in arcseconds. The systemic LSR radial velocity of the ambient molecular cloud is about 60.0 km s"' . The synthesized beam of 
the line image is 0.58" X 0.43" with a P.A. of 54.8°. The spectral resolution was smoothed to 1 km s"'. 
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W51 North 




Water Masers 



Fig. 5. Artist's conception of the molecular ring, the bipolar outflow and the central massive disk found in the high-mass protostar W51 

North. 
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